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Abstract 

Diffractive systems are quantum-mecfianical models witti point-like 
singularities where usual semiclassical approximation breaks down. 
An overview of recent investigations of such systems is presented. The 
following examples are considered in details: (i) billiards (both inte- 
grable and chaotic) with small-size scatterers, (ii) pseudo-integrable 
polygonal plane billiards, and (iii) billiards with the Bohr-Aharonov 
flux lines. First, the diffractive trace formulas are discussed with par- 
ticular emphasis on models where the diffractive coefficient diverges 
in certain directions. Second, it is demonstrated that the spectral 
statistics of diffractive models are different from the statistics of both 
integrable and chaotic systems. The main part of the paper is devoted 
to analytical calculations of spectral statistics for certain diffractive 
models. 
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1 Introduction 



Quantum chaos (or chaology) is a part of theoretical physics whose aim is the 
investigation of (multi-dimensional) quantum problems in the semiclassical 
limit h ^ 0. The modern development of this field is based on the semiclas- 
sical representation of the Green function, G{x,x '), through the sum over 
classical trajectories connecting points x and x ' (see e.g. ||T6|) 



G{x,x')= ^ Atr{x,x')e:^^(J-Sci{x,x')) (1) 

cl. tr. 

where the action, Sci{x,x '), and the pre-factor, Atr{x,x ') are computed 
from pure classical mechanics. In such an approach one implicitly assumes 
that the classical limit of quantum mechanical problem does exist which is 
not always the case. For a long time (see e.g. [0) it was known that there 
are quantum systems without 'good' classical limit. 

This paper is devoted to the investigation of a particular type of such 
models, namely, diffractive systems whose characteristic property is the pres- 
ence of point-like singularities (real or effective) in the quantum Hamiltonian. 
The existence of singularities in classical mechanics leads to the impossibility 
of continuing classical trajectories which hit these singularities. In quantum 



mechanics the situation is less dramatic [13]- [HI- Each point-like singularity 
can be described by a diffraction coefficient, D{n,n'), which determines the 
scattering amplitude on the singularity. We normalize it in such a way that 
the semiclassical expansion of the Green function in the whole space in the 
presence of a singularity at point Xq has the following form 

G{x, x') = Gq{x, x') + — y] Go{x, {xq, n))D{n, n')Go{{xo, n'),x ') (2) 

2m r^, 

n,n' 

where Go{x,x ') is the Green function (|l]) in the absence of singularity and 
Go{x, (xo, n)) is a term in (0) corresponding to a classical trajectory starting 
from point x and ending at point xq with momentum in the direction n. 

The knowledge of the Green function permits to write down the trace 
formula for diffractive systems [^, In the simplest case, when only one 
singularity is present, the density of states can be written as a sum of three 
terms 

d{E) = d{E) + dp{E) + dd{E) (3) 
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where d{E) is the smooth part of the level density given by the Thomas- 
Fermi term (plus corrections if necessary) [|16|, dp{E) is the contribution 
from classical periodic orbits, and dd{E) is the diffractive contribution. 
For integrable billiards ^ 

A °° 1 TT IT 

dpiE) = TTzY^Yl -77=^ cos{knlp - - - -nup) . (4) 
ppon=i JknL ^ ^ 



For chaotic billiard systems |[T6[| 

/ °° 1 TT 

^^^""^ = S ^ n?i |det(M;-l)|V. - r^"^^- 

In these formulas the summation is performed over all primitive periodic 
orbits (ppo) and their repetitions, Ip is the length of the ppo, z/p is its Maslov 
index, and Mp is the monodromy matrix of the ppo. 

dd{E) in Eq. (|^) is the contribution from classical orbits (called later 
diffractive orbits) which start from the singularity and end at it. 

°° 1 5 

dd{E) = ^ ImX!<^(^i'^'i)^K' ^2)^(^2, n2)/^(fi2,n3) 

m=l fi,n' 

. ■ ■G{nm,njD{n^,ni). (6) 

Here G{n, n') is the contribution to the Green function (P from a diffractive 
orbit which starts from the singular point with momentum in direction n and 
ends at this point with momentum in direction n' . The sums in (^ can be 
transformed to 

dp{E) = -- Im Aindet(l-i^) (7) 

TT Ohj 

where the operator K is defined as follows 

{k(P){n) = ^ G{n,n')D{n',n")(f){n"). (8) 



Eqs. (13)- (iD lead to a formal 'quantization' condition of diffractive systems 
(from which energy levels can be computed) 

det(l-i:) = 0. (9) 
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The main purpose of this paper is to discuss statistical properties of energy 
levels (i.e. the spectral statistics) of diffractive systems. The plan of the pa- 
per is the following. In Section 2 integrable systems with a diffractive center 
are considered and it is demonstrated that the spectral statistics of these 
systems can be computed analytically. In Section 3 the pseudo-integrable 
polygonal billiards and in Section 4 the billiards with an Aharonov-Bohm 
flux line are investigated. These models are examples of diffractive systems 
where the diffractive coefficient diverges in certain directions which consider- 
ably complicates all calculations. First, the trace formulas for these systems 
are derived and it is demonstrated that they differ from Eq. (|^). Second, by 
combining numerical and analytical arguments we investigate their spectral 
statistics. Finally, in Section 5 it is proved that the addition of a diffractive 
center to chaotic models does not change its spectral statistics. 



2 Integrable models with diffractive center 

The simplest example of diffractive systems with a constant diffractive coeffi- 



cient has been proposed in and it consists of an integrable model (below 
we shall consider a 2-dimensional rectangular billiard) with a 5-function po- 
tential 

V{x) = X6{x-xo). (10) 
The quantization condition (|^) in this case takes a particular simple form 

AG(fo,xo) = l (11) 
or (ignoring renormalization problems 0) 

AE^^^ = l (12) 

n=l ^ 

where e„ and ipn{x) are eigenvalues and eigenfunctions of the unperturbed 
system (i.e. without the (5-function potential (0)). The natural question 
arises: what is the distribution of the new eigenvalues. En, provided the 
distributions of e„ and ipn{xo) are known. 
Let us for simplicity assume that 

(i) e„ are independent random variables (as for integrable systems) with 
a step-like common distribution 

I 0, otherwise ^^^^ 
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and (ii) |-?/'(xo)P = 1 (as for a rectangular billiard with periodic boundary 
conditions) . 

Under these assumptions Eq. (O) takes the form 



N ^ ^ 

j=i ^ - ej A 

where each term in the sum is an independent random variable. 

Let us write the exact density of solutions of this equation in the following 
form 

m-H±^yj)±^. (15) 

Representing the (5-function as a Fourier integral one can express this level 
density of through the characteristic function of Eq. ( p^ 

P^^^ = exp(za(X: -i- - h E j^A-y,. (16) 

Because all ej are assumed to be independent random variables, the correla- 
tion functions of E can, in principle, be computed straightforwardly. 
In particular, the 2-point correlation function 

i?2(^l,^2)=<p(^l)p(^2)> (17) 

can be expressed in the following form 

R2{E^,E2) = J^^^[N{f{a^,a,)r-'g{a„a2) (18) 
+ N{N - l)(/(ai, a2))^-'^i(«i, a2)^2(ai, a2)]e-'("^+"^)/\ 

where 

/(tti,a2) = / ci/i(e)exp(i— + i— ^^), 
J hii — e — e 

T / N f , / \ /• «1 . "2 N 1 

^i ai,a2 = / rf/i e exp 2- + i- )— -. 

J El — e E2 — e [Ei — ey 
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When N ^ oo the direct (but tedious) calculation of these integrals gives ^ 

/22(^i,^2) =pV2(e), (19) 

where 

r2(e) = 1 + (20) 

+ daJ da2{Jl{2^/^2 ) + Jl{2^^2 ))e-2-^("i'"2)+2*(ai+a2)^ 

JO JO 



and 



1 i 

J{ai,a2) = -(«! + 02) + 7r~^("i ~ "2) 

Z ZTTAr 



+ze*^"'+"2^[(ai - a2)G{ai, 02) - iaiJo{2y/a^ ) - y/aia2Ji{2y/aia2 )], 

roo 

G{ai, 02) = e*"2 / Jo(2v^ )e^*t^t, 



Xr pX W + E' ' ^ 21^ 

This exact formula permits, in particular, to find the limiting behavior of the 
2-point correlation function. 
When e ^ 

r2(e) Ae (21) 

where 

/■oo 7r\/3 
A = 371^ xJo((3 + (5)x) Jo^(a;)c/x = — — ^ 2.72 (22) 
JO 2 

is independent of the coupling constant A, and differs from the GOE predic- 
tion (r2(e) ^ (vrV6)e ^ 1.64e). 
When e — s> 00 

For GOE r2(e) ^ 1 - l/(7re)2. 

These calculations give rise to different generalizations. 

(i) The spectral statistics can also be computed for models where the 
residues, |'?/'(xo)P are either fixed quantity (different from 1) or independent 
random variables. In particular, for the rectangular billiard with the Dirichlet 
boundary conditions the introduction of the potential (|T^) (when coordinates 
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of the scatterer are non- commensurable with the sides) leads to the following 
behavior of the 2-point correlation function at small e 



r^ie) ^ —,\n\e). (24) 



(ii) The spectral statistics of the Bohr-Mottelson model which de- 
scribes the interaction of one level with all others can be calculated by a 
similar method. 



(iii) The asymptotic behavior of the 2-point correlation function (|2T| 



can be derived without the knowledge of the exact solution (|20| ) even for more 
general cases. 

(iv) In 1^ it was demonstrated that in the summation over periodic orbits 
of integrable systems there exists hidden saddle points whose contributions 
give correct off-diagonal terms in the product of semiclassical terms like that 
in (^). This result permits to construct a specific perturbation theory by 
which one can compute successive terms of the expansion at large e of the 
2-point correlation function for integrable systems with finite diffraction co- 
efficient by semiclassical methods. 

(v) Using the results of it is possible to obtain an analog of trace 
formulas for composite operators (build from a product of the Green functions 
([^)), e.g. for conductance fluctuations |TD[. 

3 Pseudo-integrable billiards 

An interesting and important class of diffractive systems is 2-dimensional 
polygonal billiards with angles, ctj, equal rational multiples of vr 

Oi = TT (25) 

Tli 

where rrii and rij are coprime integers. These models (called pseudo-integrable 



billiards pO|, ^) have a characteristic property that their classical trajecto- 



ries belong to a 2-dimensional surface with the genus 



where N is the least common multiplier of rzj. 
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Classical trajectories in these billiards are not defined after they hit a 
corner ( pSD with rrii ^ 1. These singular angles play the role of diffraction 
centers and the quantum diffraction coefficient can be read off from the exact 



Sommerfeld solution near a wedge with angle a [26 



Dief,9,) = - sin - 

7 7 



COS - — COS — — 

7 7 



COS COS — — 

7 7 



(27) 



where 7 = a/vr. 

An important difference between these models and the ones considered in 
the previous Section is the fact that the diffraction coefficient ( p7D diverges 
in certain directions (called optical boundaries) to compensate the discontin- 
uous behavior of classical trajectories. Consequently, diffractive orbits lying 
on optical boundaries cannot be described by Eq. (^ and, first of all, a mod- 
ification of the trace formulas is required [|ll| which can conveniently be done 



by using the Kirchoff approximation |^ valid in a vicinity of optical bound- 
aries |]lT| . The basis of this approximation is the representation of the free 
(2-dimensional) Green function, G{x, x ') as an integral over a line separating 
points X and x ' 



G(f , x')=l [G(g, X ') -^Gix, q ) - G(x, q )^G{q, x ')]dq 



d 



(28) 



where Uq is the normal to the line (parameterized by g*). 

This formula is exact provided the integration is performed over the whole 
(infinite) line. The existence of the wedge reduces the integration to a semi- 
infinite region thus producing the Kirchoff approximation to the diffraction 
problem. The detailed calculations of many particular types of diffractive 
orbits in pseudo-integrable billiards have been performed in Ref. |TT|. Here 
we present only the contribution to the density of states from diffractive 
orbits lying on the boundary of multiple repetitions of a primitive periodic 
orbit of length Iq 



dioiE) 



In 



(29) 



This contribution with fixed /q and A; ^ cxo is smaller that periodic orbit 
contribution (^ but bigger that the contribution of diffractive orbits with 
finite diffraction coefficient @. 
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The diffraction coefficient ( P7[ ) can be written near an optical boundary as 
a pole term plus a finite part, D'"^^. Eq. ( p9D corresponds to the contribution 
from the pole term. It is also possible to find analytically contributions from 
the interference of D"^^^ with the pole term 

- - TffiSC g '^J + C.C. (30) 

To derive these results it was necessary to compute analytically certain multi- 
dimensional integrals. In particular, it has been shown in that 

/•OO /-OO 

/ . . . / dxi . . . c?x„e~*^^^ = 
JO Jo 

1 poo POO tt"'/^ 

= •••/ cixi . . . rfx.e-^W = (31) 

n + lJ-oo J~oo (n + 1)'^/^ 

where $(x) = xl + (xi — ^2)^ + . . . + (x„_i — x„)^ + x^, and 

00 POO poo 

dx[ ... dy,...dyr,e''^^-'y^- (32) 
J— 00 J —00 

00 roo , (n-l)/2 n-1 i 



JO 



4 g=i Jq{n - q) 



where ^(x, y) = {x - yif + {yi - 2/2)^ + . . . + {yn-i - VnY + {Un - xf. The 
calculation of these integrals is based on the invariance of the quadratic form 
$(x) under the action of a finite group generated by reflections. 

Spectral statistics of pseudo-integrable billiards in the shape of the right 
triangle with one angle equals vr/n for different n has been investigated nu- 
merically in Ref. |^ and it was demonstrated that the energy level distri- 
bution for all n (except the integrable cases n = 3, 4, 6) differs from both, 
the Poisson distribution typical for integrable models and the random ma- 
trix distribution ||18| typical for chaotic systems. It was observed numerically 
that the spectral statistics of these billiards has the following characteristic 
properties 

1. The 2-point correlation function -R2(e) when e ^ 0, i.e. there 
exist linear level repulsion as for the Gaussian orthogonal ensemble of random 
matrices. 

2. The nearest-neighbor distribution p(s) — > e~^* when s — 00 as for the 
Poisson distribution. 
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3. The number variance — > xL, when L oo. 

4. The distribution of energy levels do not change with increasing of 
energy (up to 30000 levels). 

5. With increasing n distributions stabilize quite far from random matrix 
predictions. 

6. For = 5, 8, 10, 12 the spectral statistics is close (better than 10~^) to 
the so-called semi-Poisson distribution which has the following correlation 
functions 

R,{t) = l-e-'% p(s)=4se-^ x = \- (33) 

Even the next-to-nearest distribution for these billiards agrees quite well with 
the prediction of the semi-Poisson model 

p(2,s) = ^sV2^ (34) 

The first step of analytical investigation of spectral statistics for these 
models is the calculation of the 2-point correlation form-factor in the diagonal 
approximation P|. At small r it is necessary to take into account only the 
contribution from periodic orbits which in these cases is the same as for 
integrable billiards 

K{t) = ^-Y^^8{1-AMt). (35) 

The summation here is performed over all periodic orbits. Ip is the periodic 
orbit length, Ap is the area occupied by primitive periodic orbit family. 

The billiard in the shape of the right triangle with one angle equals vr/n 
belongs to the so-called Veech polygons for which there exists a hidden 
group structure which permits the explicit calculation of the number of pe- 
riodic orbits and surfaces occupied by them. By generalizing arguments of 



it is possible to prove that for these triangular billiards 

n + e(n) 

K{0) = — ^ 36 
3[n — 2) 

where e{n) = for n odd, e{n) = 2 for n even but not divisible by 3, and 
e{n) = 6 for n even and divisible by 3. These values of K{0) are different from 
any known distributions. The existence of non-zero K{0) leads in particular 
to the linear grows of the number variance, S^(i^) — * K{0)L when L oo. 
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The contributions from diffractive and non-diagonal terms are more dif- 
ficult to find. The main problem is the divergent character of the diffraction 
coefficient which leads to the existence of terms as in Eq, ( pOf ) which grow 
quickly with / and cannot be treated in the diagonal approximation. Such 
terms should be cancelled by other contributions and this cancellation is a 
delicate procedure. 



4 Billiards with Aharonov-Bohm flux line 

Another interesting diffractive model close to the above polygonal billiards 
is a rectangular billiard with the Aharonov-Bohm flux line This model 
is defined by the Schrodinger equation 

[E + {d, - iA,f]^ = (37) 

with the vector potential of a flux line = a/r and, say, the Dirichlet 
boundary conditions on the boundary of the billiard. 

The introduction of the flux line does not change classical trajectories. 
The only difference is that any time the trajectory encircles the flux line it 
is necessary to add the phase 27ia to the semiclassical expansion (|I]). 

The exact quantum-mechanical solution of the flux line in the whole space 
gives the value of the diffraction coefficient for the scattering on the flux line 

Bil 

mA) = ^^^/^'^-'^^^'. (38) 

cos(^Y^) 

This diffraction coefficient diverges in the forward direction and, as for polyg- 
onal billiards, the trace formula requires the careful study of multiple forward 
diffraction. This can be done, as in Section 3, by using the (generalized) Kir- 



choff approximation |TT|] . For example, the contribution to the trace formula 
from the simplest diffractive orbit parallel to the base of the rectangle (with 
2 points of forward diffraction) takes the form 



2JUa - lo) ^ 
d{E) = sin^ Tca cos{2ka) (39) 

where a is the base of the rectangle and Iq is the distance of the flux line 
from the rectangle side. Note that this expression has the same dependence 
on momentum and periodic orbit length as Eq. (pOl). 



11 



For rectangular billiard with the flux line it is also possible to compute 
the value of the 2-point correlation function at small r When the ratios 
xo/a and ?/o/^ of coordinates of the flux line to the corresponding rectangle 
sides are non-commensurable irrational numbers 

K{0) = 1 - 3a + (40) 

where a is the fractional part ofa, 0<q;<1/2 and it is symmetric with 
respect to a = 1/2 when 1/2 < a < 1. 

The numerically computed spectral statistics for this model shows con- 
siderable deviations from any known distributions. 



5 Chaotic systems with difFractive center 

Consider now a chaotic system perturbed by a point-like scatterer. New en- 
ergy levels, as before, should be computed from Eq. (|12|). The only difference 
with the case considered in Section 2 is that here we shall assume that old 
(non-perturbed) energy levels, e„, are distributed as eigenvalues of one of the 



standard random matrix ensembles 18 



P{{ek})(xY[\e,~e/ (41) 

i<j 

where /?= 1,2,4 for, respectively, GOE, GUE, and GSE cases. (We ignore 
here the one-body potential needed for the confinement of eigenvalues. One 
can e.g. assume that the eigenvalues are lying on a large radius circle.) 

Within the random matrix theory the distribution of Vk = \iP{xq)\ is 
independent of the eigenvalues and is given by 



piM) = m^y-'^'^M~Nv^)- (42) 

The knowledge of the statistical distributions of the poles and the residues 
of Eq. (|12|) permits the calculation of the distribution of new energy levels. 



The first step has been done in Ref. |^ where the joint distribution of the 
new, Ej, and the old, e^, levels has been computed 

{e.}) oc n.<.-(e.-e,0(i|;.- i^.) ,,p(_^^(^^ _ ,^)), (43) 
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where p = (3/{2\N) and, due to the positivity of |ffcp, ei < Ei < Cj+i. Here 
we assume that p > and energy levels are ordered, ei < 62 < . . . < cat- 

The resuhing distribution of the new eigenvalues is defined by the expres- 
sion 

/El fE2 pEn 

dej de2... der,P{{E^},{ek}). (44) 
-00 J El JE]si-i 



In Ref. IjT^ it was demonstrated that these integrals can be computed and 
the distribution P{{Ej}) has exactly the same form as the distribution of 
non-perturbed eigenvalues given by Eq. (|3^) 



Pi{E,})^Y[\E,-E,f. (45) 

i<j 

This result is not surprising. The random matrix theory distribution is often 
considered as the result of the action of large number of small-size scatterers. 
Consequently, the addition of one more center of scattering should not change 
the spectral statistics. 

In Ref. the contribution to the 2-point correlation form-factor from 



diffractive orbits in the diagonal approximation has been computed for f- 
dimensional chaotic systems 

^'^(^) = ^(^)'^"' / \Din,n')\'dOrtdOrt'. (46) 

OpVT^ ZTT J 

According to the above arguments this additional contribution should be re- 
moved by other (non-diagonal) terms. In Ref. |1^ it has been demonstrated 



that it is the interference between diffractive orbits in the forward direction 
and periodic orbits close to the diffraction center that exactly cancel this 
term. In the derivation of this result two main ingredients were important. 
First, the uniformity principle for periodic orbits of chaotic systems, in par- 
ticular 

? rH^*'^ - 4 / '^<'>- ^"'-'^"'-'^ (^^' 

where x{q, p) is a test function defined on a Poincare surface of section {q,p). 
(qpyPp) are coordinates of points of intersections of a periodic orbit with the 
surface of section. Tp is the periodic orbit period and S = / d^ qd^pd^E — 
H{q,p)) is the phase-space volume of the constant energy surface. 

The second important point is the optical theorem for the diffractive 
coefficient which is a consequence of the unitarity of the scattering S'-matrix 
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Im D(n, n) 



/-2 



I \D{n,n')\^dOH'. 



(48) 



Using these relations one recovers [0] the invariance of random matrix resuhs 
under the addition of short-range scatterers from the semiclassical methods. 
One can easily check that ignoring the optical theorem can change the spec- 
tral statistics. 



6 Conclusion 

In this paper a few typical examples of diffractive models has been considered. 
An integrable model with a small-size scatterer with a constant diffraction 
coefficient is the simplest and the most investigated case. Under the assump- 
tion that the unperturbed system obeys the Poisson statistics it is possible 
compute rigorously the spectral statistics of this model. The main result 
is that adding a diffractive center changes completely the spectral statistics. 
The resulting statistics is characterized by level repulsion and depends on the 
value of the diffractive coefficient. Models with finite diffraction coefficient 
permit the construction of a specific perturbation theory which allows the 
term-by term computation of expansion of 2-point correlation form factor 
K{t) (and other correlation functions as well) in powers of r. The models 
such as pseudo-integrable polygonal billiards and rectangular billiards with 
a flux line are another type of diffractive models characterized by divergence 
of the diffraction coefficient in certain directions. In this case semiclassi- 
cal contribution of diffractive orbits to the trace formula differs from the 
case of finite diffraction coefficient and the spectral statistics of these models 
characterized, first of all, by non-standard value of the 2-point correlation 
form factor at small r, K{0) < 1, and consequently the linear growths of 
the number variance. For higher order terms of the expansion of correlation 
functions, the cancellation of rapidly growing terms (as in (pQ])) is not yet 
fully understood, and calculations are in progress. A simple example of a 
such cancellation is provided by chaotic systems with a point-like scatterer. 
Though it is evident (and can rigorously be proved) that when the spectral 
statistics of unperturbed chaotic system is described by one of standard ran- 
dom matrix ensembles, the addition of a diffraction center cannot change 
the statistics, in semiclassical approach this invariance requires a compensa- 
tion between diagonal and off-diagonal terms. The important point is that 
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this strong cancellation is a general phenomenon connected mostly with the 
unitarity of quantum mechanical scattering. 

Diffractive systems are an interesting and promising class of quantum- 
mechanical models. Their properties, in general, differ from standard expec- 
tations but often are accessible to analytical calculations and they open new 
perspectives in the investigation of semiclassical limit in quantum mechanics. 

Acknowledgments 

The author is very grateful to C. Schmit, N. Pavloff, P. Leboeuf, U. 
Gerland, and O. Giraud, whose collaboration permits to obtain most of the 
results of this paper. 

References 

[1] Y. Aharonov and D. Bohm, Phys. Rev. 115, 485 (1959). 

[2] S. Albeveris et al. Solvable models in quantum mechanics (Springer, New 
York, 1988). 

[3] I.L. Aleiner and K.A. Matveev, Phys. Rev. Lett. 80, 814 (1998). 

[4] M.V. Berry and M. Tabor, J. Phys. A: Math. Gen. 10, 371 (1977). 

[5] M.V. Berry Proc. R. Soc. Lond. A 400, 229 (1985). 

[6] E. Bogomolny, U. Gerland, and C. Schmit, Singular Statistics (1997) 
unpublished. 

[7] E. Bogomolny, U. Gerland, and C. Schmit, Short-range plasma model 
for intermediate spectral statistics (1997), unpublished. 

[8] E. Bogomolny, U. Gerland, and C. Schmit, Phys. Rev. E 59, R1315 
(1999). 

[9] E. Bogomolny, Nonlinearity (2000), to be published. 

[10] E. Bogomolny, O. Giraud, and C. Schmit, Trace formulas for conduc- 
tance fluctuations (2000), unpublished. 

[11] E. Bogomolny, N. Pavloff, and C. Schmit, Phys. Rev. E (2000) to be 
published. 

[12] E. Bogomolny and O. Giraud, (2000) unpublished. 



15 



[13] E. Bogomolny, P. Leboeuf, and C. Schmit, submitted to Phys. Rev. Lett. 
(2000). 

[14] A. Bohr and B.R. Mottelson, Nuclear Structure Vol. I, 302 (W.A. Ben- 
jamin, New York, 1969), 

[15] H. Bruus and N.D. Whelan, Nonlinearity 9, 1023 (1996). 

[16] M.C. Gutzwiller, Chaos in Classical and Quantum Mechanics (Springer, 
New York, 1990). 

[17] J. B. Keller, J. Opt. See. Am. 52, 116 (1962). 

[18] M.L. Mehta, Random Matrix Theory (Springer, New York, 1990). 

[19] N. Pavloff and C. Schmit, Phys. Rev. Lett. 75, 61 (1995); 75, 3779(E) 
(1995). 

[20] P.J. Richens and M.V. Berry, Physica D 2, 485 (1981). 

[21] P. Seba, Phys. Rev. Lett. 64, 1855 (1990). 

[22] A. Shudo and Y. Shimizu, Phys. Rev. E 47, 54 (1993). 

[23] M. Sieber, N. Pavloff, and C. Schmit, Phys. Rev. E 55, 2279 (1997). 

[24] M. Sieber, Phys. Rev. E 60, 3982 (1999). 

[25] M. Sieber, J. Phys. A. 32, 7679 (1999). 

[26] A. Sommerfeld, Optics (Academic Press, New York, 1954). 

[27] G. Vattay, A. Wirzba, and P.E. Rosenqvist, Phys. Rev. Lett. 73, 2304 
(1994). 

[28] W. A. Veech, Invent. Math. 97, 553 (1989); Geom. Funt. Anal. 2, 341 
(1992). 



16 



